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The dispersion relation of the high energy optical 7r-plasmons of simple hexagonal in- 
trinsic graphite was calculated within the self-consistent-field approximation. The plas- 
mon frequency u p is determined as functions of the transferred momentum q\\ along 
the hexagonal plane in the Brillouin zone and its perpendicular component q z . These 
plasmons are isotropic within the plane in the long wavelength limit. As the in-plane 
transferred momentum is increased, the plasmon frequency strongly depends on its 
magnitude and direction ((f)). With increasing angle, the dispersion relation within the 
hexagonal plane is gradually changed from quadratic to nearly linear form. There are 
many significant differences for the 7r-plasmon dispersion relations between 2D graphene 
and 3D AA-stacked graphite. They include q»- and ^-dependence and 7r-plasmon band- 
width. This result reveals that interlayer interaction could enhance anisotropy of in- 
plane 7r-plasmons. For chosen qn , we also obtain the plasmon frequency as a function of 
q z and show that there is an upper bound on q z for plasmons to exist in graphite. Addi- 
tionally, the group velocity for plasmon propagation along the perpendicular direction 
may be positive or negative depending on the choice of qy. Consequently, the forward 
and backward propagation of 7r-plasmons in AA-stacked graphite in which the energy 
flow is respectively parallel or antiparallel to the transferred momentum, can be real- 
ized. The backward flowing resonance is an intrinsic property of AA-stacked graphite, 
arising from the energy band structure and the interlayer coupling. 
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1. Introduction 

Monolayer graphene, a two-dimensional material purely made of carbon atoms ar- 
ranged in a hexagonal lattice, was created on top of a SiC^/Si substrate. ^ Its low-energy 
electronic structure shows linear dispersion corresponding to the massless Dirac fermion. 
The density of states vanishes at Fermi energy that leads to a zero-gap semiconductor. 
Such pure 2D system are identified to display rich physical properties, e.g., the peculiar 
Landau level quantization, 2 ) and the novel half-integer quantum Hall effect. 3 ) 

When a series of parallel graphene planes are periodically stacked, the 3D graphite 
crystal is formed. Graphite is one of the most important layered systems and has at- 
tracted many theoretical 4 " 26 ^ and experimental 27-34 ) studies. In terms of stacking types, 
there are four kinds of layered graphites. They include the simple hexagonal graphite 
(AA-stacked graphite), 12 > 15 > 16 ) the Bernal graphite (AB-stacked graphite), 12 > 16 ' 20 ) the 
rhombohedral graphite (ABC-stacked graphite), 13 ' 17 ' 21 - 1 and the turbostratic graphite 
(without the periodic stacking sequence). ^ The AB-stacked graphite is the most com- 
mon and stable layered graphite. The ABC-stacked graphite naturally exists just in 
combination with the AB-stacked graphite. It also has been synthesized in lab. 20 ' 21 ) As 
for the AA-stacking sequence, it only existed in the Li-intercalated graphite before, but 
was synthesized in experiments recently. 29 ~ 31 ) By the way, there were experiments to 
stress the influences of the stacking effect on physical properties. For examples: (1) angle- 
resolved photoemission (ARPES) was used to study electron-electron correlations, 32 ) (2) 
high-resolution ARPES was performed to study quasiparticle dynamics, 33 ^ and (3) the 
measurements of inelastic X-ray scattering on the spectra near Bragg reflections. 34 * 1 

Plasmon, its quantity measured in electron-energy-loss spectroscopy (EELS) is the 
loss function, which is given by Im[— l/e(g, u)]. This measurement provides informa- 
tion on the collective excitations of the electronic system. The loss function essentially 
describes plasmon excitations that arise from intraband and interband transitions. In 
contrast to optical spectroscopy, which probes only dipole-allowed (q = 0) vertical tran- 
sition, one gets access here to the momentum- dependence of the dielectric function, i.e., 
non- vertical transitions. Consequently, high-resolution EELS could be used to study im- 
portant quantities such as the character of excitonic excitations, the spatial extension 
of excitons, and the degree of localization of electronic intraband and interband excita- 
tions. Investigations on plasmons, except for understanding the elementary excitations, 
there have been recent studies involving nanoplasmonics of graphene. 35 ) The gapless 
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manetoplasmon was obtained by controlling the external electric fields, that could be 
used to develop plasmon-enhanced imaging and sensing techniques. 36 ^ Another possible 
utilization of plasmons is applied to integrated optical circuits. 

The 7r-plasmons in monolayer and muitilayer graphene and graphite had been mea- 
sured by EELS 37-40 - 1 or optical spectra. 41 ^ 44 ) The 7r-plasmon is identified from the most 
pronounced peak in the loss function. The experimental measurements show that the 
7r-plasmon frequency is from 5.1 eV to 12.0 eV and the a + n plasmon frequency cen- 
tered at 14.5 eV for monolayer graphene. The threshold 7r-plasmon frequency increases 
with increasing stacking layer-number. The 7r-plasmon frequency of graphite begins 
at ~ 7.1 eV. Additionally the monolayer graphere exhibits a linear dispersion of ir- 
plasmons, as opposed to a parabolic dispersion for multilayer graphene. However, there 
was no detailed experimental investigation on 7r-plasmons invloving the stacking order 
and anisotropy of graphites. In fact, electronic excitations, except for the magnitude 
transferred momentum, also strongly depend on its direction, and the stacking order 
could further affect the anisotropy. 

In this paper, we employ the tight-binding method to calculate the n electronic band 
structures and use the self-consistent-field approach to evaluate the dielectric function 
for AA-stacked graphite. The electronic structure shows anisotropic characteristics in 
the hexagonal plane of the first Brillouin zone. The anisotropy of 7r-band is reflected in 
special features of the dielectric function and thereby the loss function. The interlayer 
interaction further enhances the anisotropy leading to the remarkably distinct dispersion 
relations of u p versus q between 2D graphene and 3D AA-stacked graphite. 

The rest of this paper is organized as follows. A brief description of the band struc- 
ture is given in §2. The results for our dielectric function, plasmon dispersion relation 
and loss function are presented in §3. §4 is devoted to concluding remarks. 

2. Tight-binding Method and Electronic Structure 

For AA-stacked graphite, the geometric structure is formed from periodically stacked 
monolayer graphene along the z-direction. Carbon atoms in each layer have the same 
projections on the hexagonal x-y plane. The C-C bond length is b = 1.42 A, and 
the interlayer distance is c = 3.35 A. Each hexagonal unit cell includes two atoms, 
and the atomic interactions % were obtained from the study of Charlier. 11,12 ) The 
tight-binding model with 2p z orbitals is employed to calculate the electronic struc- 
ture of AA-stacked graphite. The band structure has been discussed in detail in pre- 
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Fig. 1. (Color online) 7r-electronic structures of AA-stacked graphite and monolayer graphene are 
shown along symmetric directions. The hrst BZ in the inset. 



vious work. 23,45 ) Before a further discussion on the properties of the electronic ex- 
citations, we briefly review the important features of the band structure of mono- 
layer graphene and AA-stacked graphite which differ from those of AB-stacked and 
ABC-stacked graphite. This is done for completeness and for introducing our no- 
tation. The first BZ, shown in the inset of Fig. 1, includes the symmetry points 
r[0, 0, 0], M[2tt/(36), 0, 0], K[2n/(3b), 2n/(3VSb), 0], A[0, 0, tt/c], L[2tt/(36), 0, tt/c], and 
H[2n /(3b), 27t/(3a/36), vt/c]. The suffix indicates the value of k z (in units of vr/(c) here 
and henceforth. The band structure along different symmetry directions is shown in 
Fig. 1. . 

The 7r energy bands of monolayer graphene are shown along M — > T — > K — > M 
in the first BZ in Fig. 1. The conduction and valence bands (E c and E v ) are sym- 
metric about the Fermi energy (Ep = 0). Near the K point, two linear energy bands 
intersect at the Fermi level and are isotropic in (k x , k y ). At zero temperature, only in- 
terband transitions are allowed but do not exhibit low-frequency collective excitations. 
As the temperature is increased, the new allowed intraband transitions may induce 
low-frequency plasmons. The central-energy region, along M — > V (K — > M) direc- 
tion, has a local minimum (maximum) at the M point with energy E c,v = ±7o(±2.569 
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eV). Energy dispersions around the M point make it a saddle point, which induces the 
7r-plasmon about u p ~ 27 . As for AA-stacked graphite, the energy dispersions along 
M ->■ T ->• K ->■ M = 0) and L-^A->iJ->L(A; 2 = 7r/c) are similar to 
those of monolayer graphene. The Fermi energy is moved to 0.0208 eV, inducing the 
intraband excitations near the K (H) points. Due to interlayer interactions, the band 
structure is shifted upward [downward] for k z = [k z = ir/c] with respect to that of 
monolayer graphene, and the energy separation (E s = E c — E v ) is reduced [enlarged]. 
Additionally, E s at the saddle points is increased as the states change from M — > T and 
L — > A (0 = 0°), while they decrease as the states change from M — >■ K and L — >■ H 
(0 = 90° or 30°). The plasmons around K point was studied in Ref. 22. Here, we focus 
on the 7r-plasmons. Those are excited around the saddle points M — > L in terms of 
the calculated band structure. We shall compare results for AA-stacked graphite and 
graphene. In the following section, we investigate anisotropy of n— plasmons and the 
effects due to interlayer interaction on them. 

3. Plasmon Dispersion Relation 

3. 1 The dielectric function 

The 7r-electronic excitations, due to anisotropic 7r-bands, are described by the trans- 
ferred momentum /iq = (%m cos0, hq» sin0, hq z ) and the excitation energy Tvjj. Here, 
is the angle between the in-plane transferred momentum q|| and the r — > M direction. 
Since the first BZ is hexagonal, the range 0° < < 30° is sufficient to characterize 
the direction-dependent excitations. Namely, both and + 60° cases have the same 
electronic excitations spectra. 



At arbitrary temperature T, the dielectric function calculated for bulk graphite in 
the RPA is 46 ) 



where k and q are 3D wave vectors, v q = Aire 2 /q 2 is the bare Coulomb interaction, 
and e = 2.4 is the background dielectric constant for graphite. EV(k + q) and ^(k) 
are the state energies of final and initial states and h (ti) labels conduction or valence 
bands. 5 = 0.1 70 is the energy width due to various deexcitation mechanisms and 
f(E h (k)) is the Fermi-Dirac distribution function. In demonstrating the anisotropy of 
the dispersion relation of 7r-plasmons, we first consider inelastic scattering just involving 
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q along the hexagonal plane in the BZ, i.e., k z conserved and q z = 0. Then, we will 
turn to the case of finite q z . 

The electronic excitations due to the critical points in the energy-wave-vector space 
would cause singular behavior in e. They include the K point and the saddle point 
(near M). The imaginary (real) part e 2 (ei), evaluated at the K (M) point, exhibits 
a square-root divergence (logarithmic divergence). The M point is related to the tt- 
electronic excitations, while the K point is mainly associated with the low-frequency 
ones. When 5 approaches zero, the dip-like structure in e\ becomes a discontinuity which 
accompanies zeros of e\. In the long wavelength limit (q — > 0), the dip-like structure in 
ei, due to vertical transition, occurs around oj ~ 2 70 (not shown here). With increasing 
S and q, such a structure would become broadened and occurs at higher frequency. The 
vanishing e±, if at which e 2 is small, is associated the 7r-plasmon. That is, when Landau 
damping is weak, EELS will exhibit a prominent plasmon peak. 

The calculated g-dependent dielectric function of AA-stacked graphite is mainly 
considered at = 0° (\k x ,k y ,k z ) — > \k x + q\\, k y , k z )) and = 30° (\k x ,k y ,k z ) — > 
\k x + q\\ cos(30°), k y + q\\ sin(30°), k z )). At = 0° and q\\ = 0.2 (in units of A -1 here 
and hereafter), the real part e\ (the red solid curve) and the imaginary part e 2 (the red 
dashed curve) are shown in Fig. 2(a). For e 2 , single-particle excitations (SPE) exhibit 
two peaks at uo\ = 0.436 70 and u 2 = 2.201 70, respectively. The first peak is due 
to intraband and interband transitions near the K and H points, while the second peak 
(0*2) only comes from interband transitions closely related to the M and L points. While 
the frequency is in the vicinity of uui, e\ is changed from positive value into negative 
value that exhibits a dip-like structure. As u is increased, the value of e\ gradually 
becomes positive and has one zero at which e 2 vanishes. The zero value of t\ and the 
small €2 lead to collective plasmon excitations, as shown below. As q\\ is increased, SPE 
of 62 is mainly induced by interband excitations extended from the neighborhood of K 
(H) point to M (L) point in the BZ. Figures 2(b) and 2(c) show that e 2 , except for the 
peak structure, also possesses a shoulder structure due to broadening. Additionally, the 
frequencies corresponding to the peak positions and the vanishing of e 2 increase as q\\ 
is increased. Simultaneously, the dip-like structure of e\ is less significant and occurs at 
higher frequency as q\\ is increased. This result implies that the frequency and intensity 
of 7r-plasmons will have a strong dependence on the magnitude of q\\. 

Apart from its dependence on the magnitude of the transferred momentum, the di- 
electric function also depends on the angle 0. Figure 3 shows the g| [-dependent dielectric 
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Fig. 2. (Color online) The real (ei) and imaginary (€2) parts of the dielectric function of AA-stacked 
graphite for q z = 0, = 0° are plotted as functions of frequency for fixed in-plane transferred mo- 
mentum (a) <7|| = 0.2, (b) <7|| = 0.6, and = 1.0. For comparison, the dielectric function of graphene 
is also plotted. At gii = 0.6, both of them have the same e. 
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Fig. 3. (Color online) The same plot as Fig. 2, except that <f> — 30°. 
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functions at = 30°. It is found that e has the same special structures for both = 0° 
and = 30° at small q\\, e.g., q\\ = 0.2. However, at larger q\\, peak structures (0 = 0°) 
in e are changed into nearly shoulder structures (0 = 30°). The SPE energies at = 0° 
are larger than those at = 30°, since the former (the later) are associated with E s 
along M — y T and L — y A (M — > K and L — y H). Furthermore, e\ shows more obvious 
dip-like structure and its zeros occur at smaller frequency. These results suggest that 
the spectra of plasmons have strong direction-dependence for larger transfer momen- 
tum. Furthermore, our results fully reflect the anisotropy of band structure for wave 
vectors far away from the saddle (M) and the corner (K) points. 

In order to investigate the effects due to interlayer interaction on e, Figs. 2 and 3 
(blue curves) also present the (^-dependent dielectric function for monolayer graphene 
when = 0° and = 30°. Our results show that for small q\\,e has a weak w-dependence 
on qu . However, as q\\ is increased, some significant structures in e are generated, leading 
to the same e at q = 0.6 for both graphene and AA-stacked graphite. Furthermore, at 
larger q\\, e.g., qu = 1.0, e for graphene has more structure than that of AA-stacked 
graphite. Form these results, it could be predicted that the frequency and the intensity 
of electronic excitations due to interlayer interactions will present very different in- 
dependence between graphene and AA-stacked graphite. 

3.2 The transferred-momentum-dependent loss spectra 

The loss function, defined as Im[— l/e(q, a;)], is an indicator of the intensity of elec- 
tronic excitations and is closely related to experimental probes such as light and in- 
elastic electron scattering spectroscopies. The function Im[— 1/e] is plotted in Fig. 4(a) 
for = 0°, q z = and 5 = 0.1 70 for several values of q\\. The peak above 2.8 70 
corresponds to the interband 7r-plasmon. The shoulder structure occurring at lower fre- 
quency originates from the n — y ir* single-particle excitations. The plasmon frequency 
increases with the transferred momentum, while the intensity of the plasmon peak is 
decreased as q\\ is increased. e\ approaches zero more gradually and the singular struc- 
tures in e become less pronounced (shown in Fig. 2(a)) is the main reason, since the 
derivative of e\ versus oj is inversely proportional to the strength of collective electronic 
excitations. The loss spectra as well as dielectric functions are very sensitive to changes 
in the direction of the transferred momentum. Plasmons at = 30° are distinctively 
different from those at = 0°, the plasmon frequency increases more slowly with qn and 
its intensity seems to be g||-independent, as shown in Fig. 4(b). Compared with the 3D 
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AA-stacked graphite, 2D graphene exhibits very different (^-dependent loss functions, 
as shown in Fig. 4(c) for = 0°. With increasing gy, both the plasmon frequency and 
its intensity for graphene increase more rapidly than those of AA-stacked graphite. In 
short, the high anisotropy of the electronic structure in the hexagonal plane of the BZ 
and the interlayer interactions are directly reflected in the dielectric functions and the 
frequency and intensity of the plasmon spectra. 

The q\\ and g 2 -dependent behavior of the 7r-plasmon is important in understanding 
the 7r-band features in graphite. The strong dispersion relation of the plasmon frequency 
with q means that the plasma oscillation behaves as a propagating wave with wavelength 
2ir/q and group velocity V q ui p {q). The dispersion relation of the 7r-plasmon frequency 
with q\\ is shown in Fig. 5(a) for q z = and different angles 0. The plasma frequency 
Up (~ 2.9 70 = 7.4 eV) is finite when q\\ — > and is within the region of optical 
scattering spectroscopies. Therefore, 7r-plasmon is an optical plasmon. Our results show 
that the plasmon frequency exhibits a strong 511 dependence, and that the (^-dependence 
of Up is gradually reduced with increasing 0. This result directly reflects the 7r-band 
characteristics. We further fit the dispersion relations with the quadratic equation u p = 
B Q + Biq + B 2 q 2 . B a ~ 2.92 corresponds to the threshold 7r-plasmon frequency for 
q\\ — > 0. By varying form = 0° to = 30°, B 2 decreases form B 2 = 0.793 to 
B 2 = 0.348, while B 1 increases form B x = -0.043 to B x = 0.181. The fitting result 
reveals that the dispersion relation at small presents a quadratic feature. However, 
the quadratic form is gradually reduced into the nearly linear form with increasing 0. 
The anisotropic band structures in Fig. 1 could explain the above result, since = 0° 
and = 30° respectively represent q\\ along r — >■ M and r — > K 

Figure 5(b) reveals that the 7r-plasmons of graphene are also optical modes. 47 ) How- 
ever, their oj p -q\\ dispersion relation differs from that for the 7r-plasmons of AA-stacked 
graphite. For example, for small momentum transfer, the 7r-plasmon frequency ap- 
proaches Up ~ 2.0 70 (=5.1 eV). The threshold frequency is close to the maximum 
single-particle excitation frequency from the M point and much smaller than that for 
AA-stacked graphite. This means that the depolarization shift in graphite is larger than 
graphene because the screening of the Coulomb interaction is larger in the 3D system 
compared to that in 2D. Our calculations show that for < q\\ < 1.2, the width of 
the plasmon excitation region for = 0° is Au p ~ 1.09 70 for AA-stacked graphite, 
but it is Au p PS 2.25 70 for graphene. That is, the plasmon dispersion relation for 
graphene displays stronger g-dependence than that of AA-stacked graphite. However, 
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Fig. 4. (Color online) The energy loss spectra of AA-stacked graphite for several chosen gy at (a) 
= 0° and (b) <j> = 30°. The loss spectra of graphene for different gn at <j> — 0° are plotted in panel 
(c). 
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Fig. 5. (Color online) The qu -dependent 7r-plasmon frequency is presented for several chosen angles 
4> for (a) AA-stacked graphite and (b) monolayer graphene. 
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it has also been found that as q is increased, oo p significantly shows a stronger depen- 
dence for AA-stacked graphite than for monolayer graphene. Specifically, at q\\ = 1.2 
the difference in plasmon frequency between = 0° and = 30° is Acu p ps 0.409 7 
for AA-stacked graphite, but it is Au p ~ 0.223 70 for graphene. Additionally, our curve 
fitting results for various values of show strong linearity for the dispersion relations 
of graphene. We deduce from these results that interlayer interactions significantly en- 
hance the anisotropy of in-plane 7r-plasmons, which is unusual for optical excitations. 

In Fig. 6, the loss function for ix plasmons in AA-stacked graphite is plotted versus 
frequency for = 0° and chosen in-plane transferred momentum (a) q\\ = 0.2, (b) 
5|| = 0.6 and (c) q\\ = 1.0 and 5 = 0.1 70 for several values of q z . As q z is increased, 
the sharpness of the resonance peak is reduced until there is no plasmon observable 
plasmon resonance. All the peaks appear above 2.0 70 which corresponds to the single- 
particle interband transition energy near the shoulder for the 7r — > n* transitions. The 
plasmon frequency increases with the transferred momentum q z for fixed large q», but 
decreases with increasing q z for small q» . Additionally, the intensity of the plasmon peak 
is decreased as q z is increased. Clearly, the loss spectra are very sensitive to changes in 
the magnitude of the perpendicular transferred momentum. Plasmons at large q\\ are 
crucially different from those at small q\\, since their group velocities would have different 
signs. These results show that compared with the q\\ dependence, the q z dependence of 
the loss function is very much different for AA-stacked graphite. Furthermore, while 
there is a lower bound for the in-plane transferred momentum, there is an upper bound 
for the perpendicular component of the transferred momentum for ir plasmons, further 
illustrating their anisotropy. 

In Fig. 7, the 7r-plasmon dispersion relation of AA-stacked graphite is presented as 
a function of q z . For sufficiently large q z , e.g., q z ~ 0.53, this 7r-band induced plasmon 
barely exists when the in-plane transferred momentum q\\ = 1.0. In this case, the plas- 
mons are Landau damped by the single-particle excitations. For small q z , the group 
velocity for any of the 7r-plasmon branches depicted in Fig. 7 is very small. However, its 
magnitude generally increases, but can be positive or negative depending on the value 
of The anisotropy of the band structure in the hexagonal plane of the BZ and the 
interlayer interactions are again reflected in the dielectric functions and the frequency 
and intensity of the plasmon spectra. 

Both the stacking sequence and the dimensionality play an important role in the 
electronic excitations. The high-energy plasmons are reported for the 2D single-layer 
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Fig. 6. (Color online) The energy loss spectra of AA-stacked graphite for <f> = 0° and chosen in-plane 
transferred momentum (a) qu = 0.2, (b) q\\ = 0.6 and (c) <7j| = 1.0. As q z is increased, the sharpness 
of the resonance peak is diminished until there is no plasmon resonance, indicating a cut-off, for fixed 
<7|| , as shown in Fig. 7. 
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Fig. 7. (Color online) The dispersion relation for 7r-plasmons of AA-stacked graphite as a function 
of q z for <j) = 0° and fixed gii. Each plasmon branch has an upper cut-off q z , above which there is no 7r 
plasmon. The plasmons cease to exist beyond a cut-off q z which depends on the chosen q.i . 



graphene and multilayer graphene with the AB and ABC stacking in a recent study 
by Yuan et al. by means of similar calculated theories. 48 ^ The unique features in the 
excitation spectra due to the stacking sequence are mainly from the excitations of the 
states near the M points, but they exclude the results from the Dirac Cone Approxima- 
tion. Furthermore, the more graphene layers there are the more complicated electronic 
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structures become, and it is the same for plasmons. Comparison with 3D graphite, the 
q z dependence of plasmons is absent for 2D systems. However, the dependence on the 
number of layer exhibits the interesting and rich results. They are expected to exist 
in the AA-stacked multilayer graphene. In a nutshell, the dependences of plasmons on 
the number of layer, q z , the stacking sequence, and the dimensionality provide clear 
forms of the electronic excitations needed for the experiments and are also conducive 
to further understanding graphene systems. 

4. Concluding Remarks 

In summary, in this paper, we used the tight-binding method to calculate the 7r elec- 
tron band structure of AA-stacked graphite. The band structure is highly anisotropic 
in the hexagonal plane of the BZ. But, it has a weak dependence on k z and has a 
narrower band width along the k z direction. The dielectric function is evaluated in the 
self-consistent-field approach. The 7r-plasmons are optical plasmons and show isotropic 
behavior for small transferred momentum. As the transferred momentum is increased, 
the anisotropy of the 7r-band is transformed into features appearing in e, including a 
dip-like structure in e\ and a peak structure in e 2 . The special structures are signifi- 
cantly altered by changing the magnitude and direction of the transferred momentum. 
Therefore, the plasmon frequency and the peak position of the loss functions also show 
strong qn- and 0-dependence. There exist many significant differences for the 7r-plasmon 
dispersion relations between 2D graphene and 3D AA-stacked graphite. For the former, 
the plasmon frequency has a stronger (^-dependence, resulting in a wider 7r-plasmon 
band. However, ^-dependence of 7r-plasmons, for graphite, becomes stronger with in- 
creasing q\\. These results fully reveal the role of interlayer interaction to enhance the 
anisotropy of 7r-plasmons. We have also presented the plasmon dispersion relation as a 
function of q z for fixed q\\. We show that there is an upper bound on q z for plasmons 
to exist in graphite. Furthermore, the group velocity for plasmon propagation along q z 
may be positive or negative depending on the value of qn, as we demonstrate in Fig. 7. 
Our derived results could be verified by measurements of electron-energy-loss spectra 
or optical spectra. 
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